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Abstract 

Charge centers in ionic crystals provide a channel for elementary interaction between electromagnetic 
radiation and the lattice. We calculate the electronic ground state energies which are needed to create a 
charge center - namely a F- and a //-center. In well agreement with common understanding the F-center 
results in being accompanied by a small lattice distortion whereas the //-center is accompanied by a very 
large lattice deformation. Opposite to the common understanding the additional positive charge in the 
charge center results rather to be localized on a F 4 complex than on a F 2 -complex. From the ground 
states of the charge centers we derive binding energies, diffusion barriers and agglomeration energies for 
M-center formation. These microscopic quantities are of fundamental interest to understand the dynamic 
processes which are initiated if the crystals interact with extreme intense deep ultra violet radiation. We 
further derive the equilibrium concentrations of charge centers in grown crystals. 
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I. INTRODUCTION 



Calciumfluoride (CaF 2 ) is a key material for building refractive optical elements in the deep ul- 
tra violet (DUV) range of the electromagnetic spectrum. Especially for structuring microchips an 
optical system in the DUV range is of extreme technical importance. The reason for this is that the 
smallest possible structure size, which can be obtained by a microlithographic process is propor- 
tional to the wavelength of the radiation used for structuring. In the last years immense effort has 
been made in order to provide large single crystals of CaF 2 with high optical quality. These mi- 
crolithographic processes are performed at wavelengths of excimer laser radiation sources. These 
are mainly 193nm (ArF excimer laser), 157nm (F 2 excimer laser) and 248nm (KrF excimer laser). 
CaF 2 is chosen as an optical material due to its cubic crystal structure with (nearly 1 - 2 ) perfect opti- 
cal isotropy, due to its chemical durability and its mechanical properties which make it applicable 
for lens fabrication. 

Recently high purity single crystals have been grown with impurity concentrations below the 
ppm level. However, even in these pure single crystals a finite probability for absorption processes 
remains. Since any absorption process allows for a deposition of energy in the crystal, any energy 
deposition will in turn lead to a change of material properties. Therefore there exists strong interest 
to further our understanding of the microscopic mechanism beyond UV absorption in fluoride 
crystals. 

In the past there have been a large number of investigations on charge centers in fluoride crys- 
tals in the early 70th (see e.g.^). One of the main driving forces at that time was to find materials 
for color center lasers and hosts which show zero phonon lines for spectral hole burning. Further 
on CaF 2 is used as a model system for ionic conductivity since the F-sub lattice melts approxi- 
mately 100K below the melting point of the crystal. Since that time the quality of the crystals 
was increased dramatically to to the needs of the semiconductor industry for refractive optical 
materials in the deep UV wavelength range. 

Even in an ideal material, where no point defects are present, electromagnetic radiation can 
be absorbed with a small but finite probability via a two photon process 5 6 . Such an absorption 
process will therefore quadratically increase with increasing radiation intensity and becomes of 
strong importance with increasing radiation intensity. Such a two photon process will separate 
charges and will create an electron-hole pair, provided the sum of the energy of the two photons 
exceeds the 11.2 eV of the excitonic two particle bound state of the T-exciton 7 . Above this en- 
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ergy a large probability for absorption is present . The charges of the electron and the hole have 
further the possibility to lower their energy by distorting the surrounding lattice. This means that 
a localization of the electron-hole pair takes place together with a lattice distortion. The resulting 
state is called an F-H-pa'v^^ or a self trapped exciton (STE) 11 . This F-H-pak can recombine 
with an optically forbidden transition, which leads to the fluorescence at 278nm with a relatively 
large lifetime at room temperature of l.l/x o 9 i 12 i 13 i 14 . 

In a real material there are a number of defects present. In this case the defect concentration 
is governed by entropy. The concentration at room temperature is in this case not necessarily 
determined by its room-temperature equilibrium value. This is due to the fact that the system falls 
with respect to defect statistics out of equilibrium already at the crystallization point where defect 
diffusion becomes strongly hindered. 

At high radiation rates the creation of F-H-pairs provides the major source towards radiation 
dependent change of optical material properties of CaF 2 . Due to thermal diffusion the charge cen- 
ters can move around in the material. Above 170 K 15 the charge centers are known to be mobile. 
This mobility increases on one site the probability of recombination of F- and //-centers. On 
the other site mobile charge centers can be trapped by point defects and impurities in the mate- 
rial which can lead to room temperature stable charge centers. These temperature dependence is 
a strong hint that activation energies for diffusive transport of charge centers are of crucial im- 
portance when understanding radiation dependent change in material properties of CaF 2 . In the 
present paper we calculate energies for charge center formation and derive the necessary micro- 
scopic parameters for charge center diffusion. 

The paper is organized as follows: In sectionHHwe phenomenologically describe the interaction 
between radiation and the fluoride crystal. In section [TV| we discuss the band structure of CaF 2 
which is needed as a starting point when investigating the localized charge centers. In the next 
section |V| we report calculation on a larger unit cell up to CaiosF 2 i6 where either one fluor ion 
is removed and an electron is localized on the vacancy or where an additional fluor ion has been 
added forming an essentially positively charged complex where the hole is localized. The lattice 
deformation is discussed in subsection IV PI and the obtained electronic structure in IV Al and IV CI 
In section IV El we calculate the equilibrium concentrations of F- and //-centers in an as grown 
crystal. 
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II. INTERACTION OF UV RADIATION WITH CAF 2 



For a qualitative understanding of the dynamic processes involving charge center formation and 
charge center motion we follow the arguments of 16 . We assume that a two photon process creates 
a pair of an F-center and a //-center which are spatially uncorrected. There is further no long 
range interaction between these pairs an assumption which is physically justified by the screening 
of the coulomb potential. 

Further all interaction between the F and H centers is neglected. That means we assume an 
uniform density for both centers. This approximation gets exact in the dilute limit. This means 
that the two particle correlation functions g KiK '(r) is assumed to have the uniform value 1. As the 
only driving force for F-H-pair recombination we assume statistical fluctuations of the charge 
center positions. 

PSKO = P$ 9**{r) « pW k £ {F, H} (1) 
and for the non-diagonal correlation function we assume: 



P%{r) = P [ $Ar) P? P { u ] (2) 

With the above approximations we can directly write down a system of rate equations for the 
creation of F and H center in the ideal material. 

^ = I 2 P ~ K *e{F,H} (3) 

The first term on the right site creates pairs of F- and H charge centers. This creation goes via a 
two photon process and is therefore proportional to the squared radiation intensity I 2 . The second 
term on the right site of eq. © describes the recombination of particles. This recombination is 
proportional to the probability of finding two particles within a distance smaller than the recombi- 
nation radius r and is proportional to the effective recombination rate K. 

K = 4irr (D F + D H ) (4) 

Here Dp and D H are the diffusion constants for F-center and //-center diffusion. The temperature 
dependence of these diffusion constants is known from 17 as: 

D K = d K e~^ ( 5) 

with E F = 0.7 eV and E H = 0.46 eV. In order to obtain microscopic information on the material 
parameters needed for charge center diffusion, we do an ab initio calculation of the charge centers 



together with the resulting lattice distortions. Among other results we also confirm a well known 
result 3 that there is a strong asymmetry between the number of F-centers and //-centers in CaF 2 . 
At the melting point of the F-sublattice, where the system falls out of equilibrium with respect to 
its defect concentration there is already a reasonable number of F-centers present but there are 
practically no //-centers (see sec. IV Fl) . 

HI. METHOD OF CALCULATION 

In the present work we use a commercially available density functional package (vasp 18 ). The 
code is a plane-wave code which uses ultrasoft preudopotentials. The generalized gradient ap- 
proximation (GGA) is used to calculate the exchange correlation energy. Through most of the 
calculation we set the plane wave cutoff for the numerical evaluation of the electron density to 500 
eV. This was due to restriction of computer time. In appendix |A| we show that this choice is in 
general sufficient. For the larger clusters with 32 Ca ions and 64 F ions in section|V|we perform a 
T-point calculation, which corresponds to a k-mesh of 0.575 A -1 . The criterion for the electronic 
convergence was set to AE/E = 1CT 6 . 

IV. BAND STRUCTURE OF CAF 2 

As a starting point the band structure of CaF 2 has been calculated. The resulting bandstructure 
is in well agreement with experiments and other ab initio calculations 19 . As expected for a pure 
DFT calculation the resulting band gap is with 7.7eV too small compared with the 1 1.8eV which 
are measured experimentally. 

A screened exchange calculation improves the too small band gap. 19 Also a similar calculation 
which uses an interpolation between DFT and Hartree-Fock approximation 24 arrives at a similar 
result. Details on the band structure together with a calculation of phonon modes which are com- 
pared to neutron scattering results can be found in2!l. In 2 - 1 also mechanical properties of CaF 2 and 
other fluorites are calculated. 
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V. STRUCTURE OF F- AND //-CENTER 



In order to calculate charge centers we set up a cluster of Ca 32 F 6 4 using periodic boundary 
conditions. Due to the large size of the elementary cell calculations are performed using one 
k-point only (T-point calculation). Removing one F atom and allowing for a relaxation of the 
structure into the energy minimum gives an F-center. Adding an F atom and allowing for structural 
relaxation is a configuration of a i7-center. For the structural relaxation, which takes account of 
the lattice deformation, the size of the elementary cell was kept fix. In order to obtain charge 
neutrality in the clusters an extra positive charge is trapped near the additional F ion and an extra 
negative charge is trapped in the F-ion vacancy. Due to the periodic boundary conditions a periodic 
array of charge centers is calculated where the distance is equal to four lattice constants. Since the 
electronic wave functions of such local defects drop of exponentially there should be a small but 
finite interaction between the charge centers. A check for finite size effects due to defect-defect 
interaction is made using larger systems of Ca 10 8F 2 i6 cluster and shown in section IVTl 

A. F~ -center 

Removing one F-atom leads to a cluster with Ca32F 63 . Since the F-atom had an electron closely 
attached to it, this electron is now localized in the F-vacancy. A defect level with an exponentially 
decaying wave function is created in the band gap, where the chemical potential moves in. This 
leads to a half filled single electron level, where the electron can be excited with a photon leading 
to an absorption in the optical transmission spectrum. Such a half filled state has a spin of S = 1/2 
which makes it possible for spin resonance measurements 22 to observe the defect. As a next step 
a relaxation of the structure is performed. As a result there is only minor change in the structure. 
The electronic band structure of the resulting defect is plotted in Fig. El The resulting change in 
electronic density compared to the undisturbed configuration is shown in Fig. |3] The additional 
charge is well localized at the defect and only weak lattice distortion is seen. The symmetry of the 
F-center is tetrahedron symmetry T d - In appendix |E] we calculate the splitting of excited energy 
levels due to the tetrahedron symmetry using arguments from group theory. 
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B. M-center 



A pair of two F-center is called an M -center. Usually entropy will be the dominant driving 
force for a homogeneous distribution of F-centers. On the other hand the overall lattice distortion 
which is needed to trap an electron as a F-center can be reduced if two F-center join forming 
an M-center. By removing a second F~-ion from our cluster with periodic boundary conditions 
we create an M-center. herby the two F-center can be oriented along different directions in the 
crystal. We investigate M-center in the 100 direction, the 110 and 111 direction. In the latter 
case we have two possibilities. One where a Ca 2+ is in between the two missing F-ions and 
the second one with an empty space in between. By calculating the difference in total energy 
between an M -center (Ca 32 F 62 ) and a Ca 32 F 6 4 cluster in comparison with twice the total energy 
of a F-center configuration we obtain an expression for the agglomeration energy of two F-center 
forming an M-center. 

&E M _ agg i = E C a 32 F 62 + E C a 32 F 64 ~ ^E Caz2 F &3 (6) 

We can think of this configuration as taking two F-center infinite distance away from each other 
and comparing this state with two F-center on nearest (or next-nearest or next-next-nearest) neigh- 
bor sites. The result is summarized in table |U The configuration where the two F-center are on 
nearest neighbor sites orientized along the 100 direction is the most preferred one. 

C. i7-center 

For a configuration which is used as a starting point for the F^ + -center, a F-ion is replaced by 
a pair of F-ions where according tc 12 the line connecting the two F-atoms points along the 110 
direction. In this way a Ca32F 65 -cluster is obtained. After a relaxation of the structure a stable 
configuration with localized defect level in the band gap is obtained. The resulting band structure 
is plotted in Fig. |4| The charge distribution is plotted in Fig. |5J The additional positive charge 
is located in the center of four F ions. In this way a F^ - complex is created. This is opposite 
to the standard literature e.g!^ were the if-center is described as a F^ complex. The result 
of the calculation is very clear therefore there is little doubt on the quality of the result. In the 
70'th when the microscopic models of the charge centers in fluoride crystals have been developed, 
there was no possibility to perform such a detailed microscopic calculation. The details of the 
geometry of the four F-ions which form the F|~ center is plotted in fig. One can further see 
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that the //-center is accompanied by a very large lattice distortion. This means that a large part of 
the energy, which is needed to create a localized but unbound electron-hole pair originates from 
lattice distortion around the //-center. Therefore the main part of the energy which is gained by 
trapping an electron-hole pair as a self-trapped exciton (STE) comes from the release of the lattice 
distortion around the //-center when a F-center is close. In addition this is the reason, why it was 
not possible in the past to observe single //-centers. Usually a crystal is cooled down from the 
melting point and falls out of thermodynamic equilibrium with respect to its defect concentration 
at a temperature which is shortly below the melting point. At that temperature a large number of 
F vacancies is present but only a very small number of interstitial F-ions. This means that the 
crystal even at the slowest cooling rates never reaches its thermodynamical equilibrium at zero 
temperature, where each F~ position is occupied by a F-ion. In fact there will be a reasonably 
large number of F - vacancies. By heating the crystal and quenching it down to liquid nitrogen 
temperatures it is possible to additively color the system by creating F-centers. This argument is 
quantified in more detail in sec. IV El A defect which is discussed in CaF 2 as well is the v^-center. 
In the Vfe-center two F -ions are shifted out of their equilibrium positions forming an F^-molecule. 
We did not find a stable configuration for such a bare center and believe that it is only stabilized in 
the vicinity of dopant ions. 

D. lattice relaxation 

In the following the energy to create an electron-hole pair in the form of an F-//-pair is calcu- 
lated, this is the total energy of an F-//-pair. The undisturbed lattice has a total energy of -17.696 
eV per CaF 2 unit. For a Ca 32 F 64 this is —17.64 32 = —564.63 eV. A F-H-pak which has infinite 
distance is the electronic and fully converged structure of a Ca 32 F 6 3 cluster plus the structure of a 
Ca 32 F 65 cluster. 

AE F _ H = 2E Ca . i2F(i4 - E Ca32F(i3 - E Ca32F65 

= 2 564.63 - (-556.01) - (-565.08)eK 

= -8.17eV (7) 

via an optically forbidden transition the F-//-pair can recombine with the relative long lifetime 
of 1.1 /is and finally relax into the undisturbed configuration. The STE fluorescence which corre- 
sponds to the F-//-pair recombination is known to be at 278 nm 13 . In Fig. |7]we show a schematic 
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plot of the energy scheme which relates the measured 278 nm which is a photon energy of 4.463 
eV to the calculated energy difference of 8.17 eV. 



E. Equilibrium concentrations of charge centers in a grown crystal 

In subsections IV Al and IV CI we saw that there is a strong asymmetry in the energy which is 
needed to produce a F- or if-center. The main reason for this asymmetry lies in the strong lattice 
distortion which results when adding a F-ion into an interstitial lattice site in order to create a H- 
center. We now estimate the concentration of F- and //-center which is obtained when the system 
falls out of equilibrium. Therefore we compare the total energy of Fluorine in the four possible 
configurations, (i) as Flour in the regular CaF 2 lattice, (ii) as Flour in an interstitial lattice site, 

(iii) the energy gain which is obtained when removing a F ion in order to create an F-center and 

(iv) the Fluor in the F 2 atom of gaseous Fluor. From section IVTJ1 we know that the energy which 
is needed to create a F-H -pair is 8.17 eV. We now ask how this energy is distributed between the 
F- and the //-center. For the F-center we obtain: 



Which is the energy gain when creating the F-center and forming half of a F 2 molecule. In** a 
much larger value («7.9 eV) is obtained for F _ -center formation by comparing to the energy of 
a bare F atom instead of a F 2 molecule. This is in agreement with older ab-initio studies 25 , where 
small clusters of only 18 ions where studied and a diffusion activation energy for the F-center of 
1.69 eV results. For the //-center we obtain in the similar way: 



The sum of both energies gives 8.023 eV which is the energy needed to create an F-H -pair. In prin- 
ciple one should add in Eq. © and subtract in Eq. © the chemical potential of a corresponding 
F 2 gas phase, which in general depends on temperature and pressure AE F — > AE F + fi(T,p, N) 
and AE H — > AE H — fi(T,p, N). For simplicity and due to a lack of knowledge of the detailed 




1.35eV 



(8) 




( Ec a32 F 6 4 + _ Eca 32 F 65 J ~ AE F _ H 

(-564.63 + (-1.79) - (-565.08)) - 8.17e^ 



6.82eV 



(9) 
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environment, we set this term to zero. The point where the system falls out of equilibrium with 
respect to its F-ion concentration is given by the melting temperature of the Flour sub-lattice in 
CaF 2 . This occurs at approximately 1500K which is 150K below the melting point of the crystal. 
At this temperature the equilibrium concentration of F-center is: 

1 AE F 

(nf ) = — e'^ = 2.5 10 23 m" 3 = 2.5 10 17 cm" 3 (10) 
and for the iJ-center also at a temperature of 1500 K it is: 

1 AE H 

( n J) = — e"^ = 1.2 10 4 m~ 3 = 0.012 cm" 3 (11) 

where v = a 3 = (5.46 10~ 10 ) 3 m 3 = 1.67 10~ 28 m 3 is the volume of the elementary cell. This 
result means that in ideally grown crystal there is already a strong asymmetry between Fluor va- 
cancies and Fluor interstitial. There are practically no i/-centers but already a finite number of 
F-centers. Using the molar weight, the density of CaF 2 and the Avogadro constant, NA, we can 
calculate the molar concentration of charge centers in thermal equilibrium at different tempera- 
tures and as well an absolute particle number in ppm. This is done in table HU To calculate the 
equilibrium densities in ppm we used 

« /F ) Ppm =« F )^f (12) 

where u Ca F 2 — 78.08 g/mol is the molar weight of CaF 2 , p = 3.31 g/cm 3 the density of CaF 2 and 
NA the Avogadro number. 

There is even a very small but finite change in the density of the material due to the equilibrium 
concentration of F-vacancies. This change in density is: 

(n F ) 

An = 19.00 g/mol 

NA 

9 ^ in 17 

= 19.00 - g/cm 3 

6.022 10 23 S 

= 7.9 10~ 6 g/cm 3 = 7.9 10~ 3 kg/m 3 (13) 

From that follows that an extremely inhomogeneous distribution of charge centers can in principle 
have an influence on the density homogeneity of the material and therefore also on the refractive 
index homogeneity. Since, however, the bare charge centers are mobile at room temperature, this 
process will lead to a homogeneous distribution already. 
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F. Finite size effects 



In order to estimate finite size effects we configured a larger cluster with CaiosF2i6 f° r the 
unperturbed configuration. In this cluster the resulting configuration for the F-center and for the 
iJ-center from the smaller clusters was inserted. For these calculations we restricted ourself to 
a plane wave cutoff of 400 eV only due to limitations of computer power. After this first one 
electronic configuration was performed and the forces on the ions where calculated. As a result 
there is still a non negligible finite size effect present. For the F-center we obtain: 

A T?F elec— relax _ _ „ rpelec— relax 

l\.Z fin _ size - lbJ^CaF 2 + ^Ca 32 F 63 ~ &Ca W8 F 215 

= -1902.59- (-1908.57) 

= 5.98eV (14) 
If we allow for a structural relaxation in the larger cluster we get: 

A rpF struct— relax rpstruct— relax 

l\.£j fin _ size - (b&CaF 2 + ^Ca S2 F &3 ~ ^Ca 10S F 215 

= -1902.59 - (-1910.42) 

= 7.83eV (15) 

Where E Ca F 2 is the total energy of a single CaF 2 unit cell and E Ca32 F e>3 is the total energy of the 
cluster from subsection IV Al and Ec ai0S F 215 is the total energy of the large superstructure. The 
differences in energy result from an interaction between a single charge center and its neighbors 
due to the periodic boundary conditions of the underlying system. The major part (5.98 eV) stems 
from the electronic interaction wheres only a small part (7.83 eV- 5.98 eV = 1 .85 eV) stems from an 
interaction between the different lattice distortions around each individual charge center. Therefore 
the major contribution of the relatively large energy difference originates from the long range 
Coulomb interaction between the charge center and it's neighbors due to the periodic boundary 
conditions. 

For the -center we obtain from a similar calculation for the electronic part only: 

A r?H elec— relax _ „ „ rpelec— relax 

l\ti fin _ size - C0ti C aF 2 + &Ca 32 F 65 ~ &Ca 10s F 21T 

= -1911.67- (-1912.82) 

= l.lbeV (16) 
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And for the structural relaxation we get: 



AF 



struct— relax 



7QE C aF 2 + Eca- i2 F 65 ~ E\ 



^struct— relax 
0(1108^217 



fin— size 



1911.67- (-1913.13) 



1.39eV 



(17) 



In this case as well the electronic part (1.15 eV) as well as the structural part due to the interaction 
between lattice distortions (1.42 eV - 1.15 eV = 0.27 eV). 

The reason why finite size effects are much smaller for the ff -center is due to the more efficient 
screening of the Coulomb charges in the case of the if -center. In this case there is a much larger 
but also much stronger localized lattice distortion which leads to more strongly localized electronic 
wavefunctions of the defect level. This can be also seen by the weak but still present dispersion of 
the defect level in Figs. EJand|U 

VI. CONCLUSION 

In the present paper we calculate the ground state of the trapped electron (bare F-center) and 
of the trapped hole (bare if -center) in CaF 2 . These charge centers are the elementary excitations 
which are needed to understand the interaction between CaF 2 and DUV radiation. Since DUV 
radiation can even in an ideal crystal excite electron-hole pairs via a two photon process, the elec- 
tron or the hole can localize, gaining energy via a lattice distortion. The formation energy for 
a F-H -pair is, however, unevenly shared between the electron and the hole. This is the reason 
why already in an as-grown-crystal a finite number of F-centers of the order of 5 ppm is present, 
while there are no if -centers present. This explains why in Fluorite crystals there are F-centers 
observed but ff -centers are observed only if they are stabilized by impurities. Calculating the lat- 
tice distortion around the F-center, we found no surprise. However when investigating the lattice 
distortion around the if -center we found that the state of smallest total energy is a configuration 
where the direction connecting the two interstitial F ions points along the 1 10 direction of the crys- 
tal. Further the charge is distributed rather around in a F^~ complex rather than in a F^ complex. 
For the important problem of radiation stability of CaF 2 under intense excimer laser radiation, the 
trapping of these charge centers on impurity ions has to be considered. The radiation stability of 
CaF 2 is of extreme importance for the optical microlithography for structuring semiconductors. 



12 



APPENDIX A: TEST OF NUMERICAL ACCURACY DUE TO PLANE WAVE CUTOFF 

In a strongly ionic crystal like CaF 2 there is a large spatial fluctuation of the electron charge. 
Therefore it is expected that a large number of plane waves (in k-space) has to be used to achieve 
a numerical accuracy which allows quantitative predictions. 

In Fig. [8] we have plotted the total energy of a CaF 2 system as a function of the plane wave 
cutoff. In Fig. |9]we show the resulting lattice constants after a structural relaxation. For cutoff 
energies of 500 eV and larger the system starts to saturate. To minimize calculation time we 
have used a plane wave cutoff of 500 eV throughout the calculations of the larger clusters for the 
calculation of charge centers and impurities in CaF 2 . In Fig. |5]the resulting lattice constant of 
CaF 2 is shown as a function of the cutoff energy. 

APPENDIX B: SPLITTING OF ENERGY LEVELS OF CHARGE CENTER DUE TO SYMME- 
TRY CONSIDERATIONS 

In this appendix we give a few thoughts about symmetry, which are based on group theoretical 
arguments. We assume that a charge center forms a kind of hydrogen atom which is perturbed 
by the symmetry of the underlying crystalline lattice. In the case of the F-center this is obvious 
since the fiuor vacancy forms a net positive charge in which the excess electron is localized due 
to Coulomb attraction. The point symmetry group of the F-center is tetrahedron symmetry, T<f. 
In the case of the if -center this is less obvious. The excess fluorine ion brings a negative charge 
which is compensated by a positive charge which is localized among the four F-ions as shown in 
fig - HD The point symmetry of the whole if-center is C 2 „. 

The optically allowed transitions within the hydrogen atom are given by well known rules 
which follow from symmetry considerations in the rotation group, 0(3)=R+Q. The angular part of 
the eigenfunctions in the 0(3) group are the spherical harmonics, Df . Where I is the angular index 
of the spherical harmonics and p the parity (+ or -) of the wave function. Under the perturbation 
of the external field these (21+1) times degenerate eigenvalues split in a well defined way. These is 
given in correlation tables which can be found in standard textbooks of group theory 26 . In table IUTI 
we give the part of the correlation table which is relevant for the symmetry of the charge centers. 
From tableHniwe learn the following important things: (i) the splitting of the energy levels in the 
excited states and second the numbers of independent parameters which are required to expand the 
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crystal field up to fourth order in I. In the case of the F-center only the second excited state splits 
into a twofold (E) and a threefold (T 2 /i) degenerate energy level. The third excited state splits 
into three levels, one highly symmetric A 2 /i and two levels which are each threefold degenerate 
Ti and T 2 . In the case of the if -center all degeneracies up to fourth order are canceled, already 
the first excited state splits into three non degenerate energy levels A 2 /i, B\ and B 2 . The second 
excited energy level splits into five different levels and so on. By (ii) counting the number of Ai 
representations we know the number of independent parameters which are needed for an expansion 
of the crystal field up to fourth order in /. For the F-center these are only 3 parameters whereas 
for the //-center with its much lower symmetry 1 1 parameters are required. 
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figure captions 

FIG. 1: band structure of CaF2 from an LDA calculation. As expected the band gap results in a too small 
value, whereas the main features like e.g. the indirect transition are obtained in well agreement with exper- 
iments. 

FIG. 2: band structure of Ca32F@3 from an LDA calculation. A weakly dispersive defect level shows up in 
the band gap. the Fermi level marked as of the energy scale is in the half filled defect level indicating a 
free spin of the defect. 

FIG. 3: Change in electronic density due to formation of an F-center. The F-center is well localized and 
accompanied by only a very small lattice distortion. 
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FIG. 4: band structure of Ca32Fg3 from an LDA calculation. A weakly dispersive defect level shows up in 
the band gap. the Fermi level marked as of the energy scale is in the half filled defect level indicating a 
free spin of the defect. 

FIG. 5: Change in electronic density due to formation of an H-center. The -center is well localized and 
accompanied by a huge lattice deformation. The additional positive charge is located at the center of four 
F~-ions forming an F^ - complex. The symmetry of the //-center is C2 V - In appendix |B] we estimate the 
level splitting due to symmetry arguments. 

tables 



type of 
M-center 


AEM-aggl 

/eV 


M-100 


-0.85 


M-110 


-0.14 


M-lll (Ca) 


-0.31 



TABLE I: Agglomeration energy for the formation of an M-center in CaF2. According to our result in the 
most favoured configuration the two F-vacancies are oriented along the 100 direction. In this direction the 
gain due to lattice relaxation is largest. 
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FIG. 6: The geometry of the F4 - , ii-center as it lies in the (-1,1,0) plane is shown. The angular as they 
result from the calculation are plotted and given in ° in fig. a). The distances between the four F-ions are 
given in A. The four F-ions form a perfect plane above and below which a Ca 2+ -ion is located. The overall 
symmetry of the if -center obtained via our calculation is C2 V as can be seen from fig. b). 

FIG. 7: Schematic plot of the energy levels of the F-ii-pair before and after recombination. Resulting from 
the calculation there is a large part of the energy stored in a lattice deformation (Stokes shift). Therefore the 
measured 278nm fluorescence is smaller than the calculated energy difference of 8.17 eV. 

FIG. 9: The lattice constant of CaF2 after structural relaxation is plotted as a function of plane wave cutoff. 



dopant 


^Edefect 

/eV 


n eq (T ~- 
/m~ 3 


= 300^) 
/ppm 


n eq (T = 
/m~ 3 


= 1500K) 
/ppm 


n eq (T = 
/m~ 3 


= 1650K) 
/ppm 


F-center 


-1.35 


2.49 10 4 


9.74 i<r 19 


1.29 10 23 


5.04 


3.43 10 23 


13.4 


if-center 


-6.82 


1.75 1(T 91 


6.84 i<r 114 


1.20 10 4 


4.70 10~ 19 


1.72 10 6 


6.72 10~ 17 



TABLE II: For F- and if-center the equilibrium concentrations at different temperatures are calculated. 
The results are given as well in absolute concentrations [cm -3 ] as well as in molar fractions [ppm]. 



FIG. 8: The total energy for CaF2 as a function of plane wave cutoff. For cutoff energies of 500 eV the 
system starts to saturate. The upper (red) curve shows the result without structural relaxation and the lower 
(green) curve shows the result with structural relaxation of the elementary cell. 
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F-center 


i^-center 


0(3)=R x Q 






D ± 
u 


A1/2 


A1/2 


Df 

1 




Ao / 1 +Bi +B2 




E+T 2 /i 


A 1 / 2 +2Bi+2B 2 




A 2/1 +T 1+ T 2 


A 1/2 +2A 2/ i+2B 1 +2B 2 




A 1/2 +E+Ti+T 2 


3A 1/2 +2A 2/1 +2Bi+2B2 









TABLE III: For the first four excitations of a charge center we give the relevant part of the correlation table 
for a splitting of the energy levels with respect to the symmetry of the F- and the H- center. 
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